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Abstract 

In this paper, we study the possibihty of building Yang-Mills(YM) field dark energy models with 

OO . equation of state (EoS) crossing -1, and find that it can not be realized by the single YM field 

O ■ 

' models, no matter what kind of lagrangian or initial condition. But the states of —1 < u; < and 

(N ; 

J> ' u; < — 1 all can be naturally got in this kind of models. The former is like a quintessence field, and 

o : 

■ the latter is like a phantom field. This makes that one can build a model with two YM fields, in 

\ which one with the initial state of — 1 < lj < 0, and the other with cj < — 1. We give an example 

\ model of this kind, and find that its EoS is larger than -1 in the past and less than -1 at the present 

\^ • time. We also find that this change must be from a;>— lto<— 1, and it will go to the critical 

in '. 

^T) . state of cj = — 1 with the expansion of the Universe, which character is same with the single YM 

o : 

T— I ' field models, and the Big Rip is naturally avoided. 

in : 
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I. INTRODUCTION 

Recent observations on the Type la Supernova (SNIa)|l|, Cosmic Microwave Background 
Radiation (CMB)|^ and Large Scale Structure (LSS) sl all suggest that the Universe mainly 
consists of dark energy (73%), dark matter (23%) and baryon matter (4%). How to under- 
stand the physics of the dark energy is an important mission in the modern cosmology, which 
has the EoS of < —1/3, and leads to the recent accelerating expansion of the Universe. 
Several scenarios have been put forward as a possible explanation of it. A positive cosmo- 
logical constant is the simplest candidate, but it needs the extreme fine tuning to account 
for the observed accelerating expansion of the Universe. This fact has led to models where 
the dark energy component varies with time, such as quintessence modelsyi, which assume 
the dark energy is made of a single (light) scalar field. Despite some pleasing features, these 
models are not entirely satisfactory, since in order to achieve Qde ~ (where Qde and Qm 
are the dark energy and matter energy densities at present, respectively) some fine tuning 
is also required. Many other possibilities have been considered for the origin of this dark 
energy component such as a scalar field with a non-standard kinetic term and k-essence 
modelsj^, it is also possible to construct models which have the EoS oi uo = p/ p < — 1, 
the so-called phantomlsl. Some other models such as the generalized Chaplygin gas (GCG) 
models the vector field models^ also have been studied by a lot of authors. Although 
these models achieve some success, some problems also exist. One essential to understand 
the nature of the dark energy is to detect the value and evolution of its EoS. The observation 
data shows that the cosmological constant is a good candidate jol, which has the effective 
equation p = —p, i.e. uj = —1. However, there is an evidence to show that the dark energy 
might evolve from > — 1 in the past to < —1 today, and cross the critical state of 

= — 1 in the intermediate redshift [10| . If such a result holds on with accumulation of 
observational data, this would be a great challenge to the current models of dark energy. It 
is obvious that the cosmological constant as a candidate will be excluded, and dark energy 
must be dynamical. But the normal models such as the quintessence fields, only can give 
the state of — 1 < c<j < 0. Although the k-essence models and the phantom models can get 
the state of u < —1, but the behavior of u crossing —1 can not be realized, and all these 



will lead to theoretical problem in field theory. To answer this crossing phenomenon o 



lot of people have advised some more complex models, such as the quintom models 11 
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which is made of a quintessence field and a phantom field. The model with higher derivative 
term has been suggested in Ref. 3], which also can get from co'>— ltoco'<— 1, but it also 
will lead to theoretical difficulty in field theory. 

We have advised that the YM field|l4, 15] can be used to describe the dark energy. There 
are two major reason that prompt us to study this system. First the normal scalar models 
the connection of field to particle physics models has not been clear so far. The second 
reason is that the weak energy condition can not be violated by the field. The YM field 
we have advised has the desired interesting featured: the YM field are the indispenable 
cornerstone to any particle physics model with interactions mediated by gauge bosons, so 
it can be incorporated into a sensible unified theory of particle physics. Besides, the EoS of 
matter for the effective YM condensate is different from that of ordinary matter as well as 
the scalar fields, and the state of — 1 < u; < and uj < —1 can also be naturally realized. 
But if it is possible to build a YM field model with EoS crossing —1? In this paper, we 
focus on this topic. First we consider the YM field with a general lagrangian, and find the 
state of ~ — 1 is easily realized, as long as it satisfies some constraint. From the kinetic 
equation of the YM field, we find that u +1 oc a^^ with the expansion of the Universe. But 
no matter what kind of lagrangian and initial condition we choose, this model can not get 
a behavior of u crossing —1. But it can be easily got in the models with two YM fields, one 
with the initial condition of > —1, which is like a quintessence field, and the other with 
UJ < —1 like a phantom field. 

This paper is organized as follows. In section 2 we discuss the general YM field model, 
and study the evolution of its EoS by solving its kinetic equation. But we find that this kind 
of model can not get the state of u crossing —1. Then we study the two YM fields model in 
section 3, and solve the evolution of uj with scale factor for an example model. We find that 
u crossing —1 can be easily realized in this model, which is very like the quintom models. 
At last, we have a conclusion and discussion in section 4. 



II. SINGLE YM FIELD MODEL 



In the Ref. 151], we have discussed the EoS of the YM field dark energy models, which 
has the effective lagrangian ly, [17 ] 



'eff 



F/2g\ 



(1) 



here F = —{1/2)F^^F°'^^ plays the role of the order parameter of the YM condensate, and 
g is the running coupling constant which, up to 1-loop order, is given by 

1 F 

Mn|--l|. (2) 

Thus the effective lagrangian is 

Le// = ^i^ln||^l, (3) 
where e ~ 2.72. b = llA^/247r^ for the generic gauge group SU{N) is the Callan-Symanzik 



coefficient 



18|, K is the renormalization scale with the dimension of squared mass, the only 



model parameter. The attractive features of this effective YM action model include the 
gauge invariance, the Lorentz invariance, the correct trace anomaly, and the asymptotic 
freedom With the logarithmic dependence on the field strength, i-ieff has a form similar 
to he Coleman- Weinberg scalar effective potentialjl^, and the Parker- Raval effective gravity 
lagrangian {20!. 

It is straightforward to extend this model to the expanding Robertson- Walker (R-W) 
spacetime. For simplicity we will work in a spatially flat R-W spacetime with a metric 

ds"^ = a^{T){dT'^ - 'jijdx'dx^), (4) 

where we have set the speed of light c = 1, ■jij = 5* denoting the background space is flat, 
and T = J {ao/a)dt is the conformal time. Consider the dominant YM condensate minimally 
coupled to the general relativity with the effective action, 

S = I ^[-Y^ + Le//]rfV (5) 

where g is the determinant of the metric (7^1^. By variation of S with respect to the metric 
gf^", one obtains the Einstein equation G^i, = SvrGT^i,, where the energy-momentum tensor 
is given by 

^/^^ = E ^2F:sF'"'' + eF;,Fr. (6) 

The dielectric constant is defined by e = 2dLeff /dF, and in this one- loop order it is given 
by 

6 = Mn|^|. (7) 

This energy-momentum tensor is the sum of the several different energy-momentum tensors 
of the vectors, T^,y = Yla ^"'''F^^, neither of which is of prefect-fluid form, which can make the 



YM field being anisotropy. This is one of tlie most important cliaracter of the vector field 
dark energy models j^. If it is true and this anisotropy YM field is dominant in the Universe, 
this will make the Universe being anisotropy, one would expect an anisotropy expansion of 
the Universe, in conflict with the significant isotropy of the CMBj^l]. But on the other hand 
there also appear to be hints of statistical anisotropy in the CMB perturbations 22] . But 
here we only consider the other case. For keeping the total energy-momentum tensor T^,^ is 
homogeneous and isotropic, here we assume the gauge fields are the functions of only time t, 
and = |craAj^(t) (here cTq are the Pauli's matrices) are given by = and = 61A(t). 
Define the YM field tensor as usual: 



f;, = d,Al - d^Al + r'^A\Al 



(8) 



where / is the structure constant of gauge group and / = e for the SU (2) case. This 
tensor can be written in the form with the electric and magnetic field as 

/ Ex E2 



V 



(9) 



— E\ -B3 — B2 
—E2 — -B3 Bi 
\-Es B2 -Bi / 

It can be easily found that Ef = E2 = E^, and Bf = B2 = B^. Thus E has a simple form 
with F = E"^ — B^, where E"^ = Yl^=i ^1 ^'^'^ ~ S^=i ^I- ■'■^ ^^^^ case, each component 
of the energy-momentum tensor is 



^^X = ^^{B^-E^)6l+'-EX 



A" 



{a)rpi 
3 



^{B^ - E^)5] + ^-EH]5'; - \B^5i{l - 6]) 



(10) 



(11) 



Qg^' ' ' 3 ' ' 3 

Although this tensor is not isotropic, its value along the j = a direction is different from 

the one along the directions perpendicular to it. Nevertheless, the total energy-momentum 

tensor T^^ = X]a=i ^"'^T^u has isotropic stresses, and the correspondiiig energy density and 

pressure are given by (here we only consider the condition oi B^ = 0) 15] 

E\ E^ - 



and its EoS is 



UJ 



e-36 
3e + 36' 



(12) 
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It is easily found that at the critical point of e = 0, which follows that u = —1, the Universe 
is exact a de Sitter expansion. Near this point, if e < 0, we have u; < — 1, and e > follows 
uj > —1. So in these models, the EoS of < u; < — 1 and u; < — 1 all can be naturally 
realized. 

For studying the evolution of this EoS, we should solve the YM field equations, which is 
equivalent with solving the Einstein equation[l5|. By variation of S with respect to A^, one 
obtains the effective YM equations 

d^ia^e F^^") + r^^A^^ia^e F^^"") = 0. (14) 

For we have assumed the YM condensate is homogeneous and isotropic, from the definition 
of -F^^, it is easily found that the u = component of YM equations is an identity and the 
i = 1,2,3 spatial components are: 

dria^eE) = 0. (15) 
If e = 0, this equation is also an identity. When e 7^ 0, this equation follows that [fsj] . 

(3 0, a-', (16) 

where we have defined /3 = e/b, and used the expression of e in Eq.Q. In this equation, 
the proportion factor can be fixed by the initial condition. This is the main equation, which 
determines the evolution of f3, and /3 directly relate to the EoS of the YM field. Combining 
the Eqs. fll3l) and (ITGll . one can obtains the evolution of EoS in the YM field dark energy 
Universe. In Fig. [1] , we plot the the evolution of u in the YM field dark energy models with 
the present value ujq = —1.2 and tuo = —0.8, and find that the former one is very like the 
evolution of the phantom field, and the latter is like a quintessence field. They all have same 
attractor solution with uj = —1. So in these models, the Big Rip is naturally avoided. This 
is the most attractive feature of the YM field models. 

In the Eq.( [T6i) . the undetermined factor can be fixed by the present value of EoS uq, 
which must be determined by observations on SNIa, CMB or LSS. In this paper, we will 
only show that the observation of CMB power spectrum is an effective way to determine it. 
The dark energy can influence the CMB temperature anisotropy power spectrum (especially 
at the large scale) by the integral Sachs- Wolfe (IS W) effect [23;] • Consider the fiat R-W metric 
with the scalar perturbation in the conformal Newtonian gauge, 

ds^ = a\T)[{l + 2(P)dT^ - (1 - 2i!)-i,jdx'dx^]. (17) 



The gauge-invariant metric perturbation ifj is the Newtonian potential and is the pertur- 
bation to the intrinsic spatial curvature. Always the background matters in the Universe 
are perfect fluids without anisotropic stress, which follows that (p = ip. So there is only one 
perturbation function in the metric of f[T71) . and its evolution is determined by [24] 



f + m{l + - ^VV + [(1 + 3^)H' + 27i']<P = ATxGa\5p - ^(5p), (18) 



where Ti = a'/a, and the 'prime' denotes d/dr. The pressure p = YliPi^ energy density 
p = YliPi^ which should include the contribution of baryon, photon, neutron, cold dark 
matter, and the dark energy. Especially at late time of the Universe, the effect of the dark 
energy is very important. We remind that the ISW effect stems from the time variation of 
the metric perturbations, 

/^u rxLss 
y[y^ dxi^' + ^%ikx)r, (19) 

where Xlss is the conformal distance to the last scattering surface and ji the I'th spherical 
Bessel function. The ISW effect occurs because photons can gain energy as they travel 
through time-varying grav itational wells. One always solves the CMB power spectrum in 

261]. In Fig. [2], we plot the CMB power spectrum at large scale 



the numerical methods 



25 



with these two kind of YM dark energy models, where we have chosen the cosmological 
parameters as: the Hubble parameter h = 0.72, the energy density of baryon Q^h'^ = 0.024, 
and dark matter Qdmh"^ = 0.14, the reionization optical depth r = 0.17, the spectrum index 
and amplitude of the primordial perturbation spectrum being Ug = 0.99 without running and 
A = 0.9. Where we haven't consider the perturbation of the dark energy. From this figure, 
one can find that the values of the CMB power spectrums are very sensitively dependent 
on uq. Comparing with the ACDM model (which is equivalent with the YM model with 
uq = —1), the model with t^o = —0.8 > —1.0, which is like the quintessence field model, 
the CMB spectrums have smaller values, especially at scale of / < 10, the difference is very 
obvious; but the model with c^o = —1.2 < —1.0, which is like the phantom field model, the 
CMB spectrums have larger values. For the evolution of EoS is only determined by the uq, 
the value of it can be determined by fitting the CMB observation. It is obvious that the 
recent observations on the CMB power spectrums at large scale from WMAP satellite have 
large error. The further results will depend on the observation of the following WMAP and 
Planck satellites. 
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Now let's return to the evolution of uj. From Fig.[l], one finds that uj crossing —1 can 
not be realized in these models with a single YM field, no matter what values of we have 
chosen. For studying it more clear, assume the YM field has an initial state of |u; + 1| -C 1, 
which follows that /3 ^ 1, the Eq. flTB]) becomes 

(3 oc (20) 

The value of j3 will go to zero with the expansion of the Universe. This means that E will 
go to a critical state of E"^ = k^. And the EoS is 

+ (21) 

3 

This result has two important characters: i) with the expansion of the Universe, uj will go 
to the critical point oi uo = —1. This is the most important character of this dark energy 
model, which is very like the behavior of the vacuum energy with u; = — 1; ii) the value of 
uo > —1 and < — 1 all can realized, but it can not cross —1 from one area to another. This 
character is same with the scalar field such as the quintessence field, the k-essence and the 
phantom field models. 

It is interesting to ask: if these characters are correct just for the YM model with the 
lagrangian as formula ([3])? Whether or not one can build a model, whose EoS can cross —1? 
So let's consider the YM field model with a general effective lagrangian as: 

= G(F)F/2, (22) 

where G{F) is the running coupling constant, which is a general function of F. If we choose 
G{F) = 6 In 1^1, this effective lagrangian returns to the from in Eq.([3]). The dielectric 
constant also can be defined by e = 2dLeff/dF, which is 

e = G + FGp. (23) 

Here Gp represents dG/dF. We also discuss the homogeneous and isotropy YM field with 
electric field (5 = 0), then the energy density and the pressure of the YM field are: 

P = E\e-^), (24) 
P=-E\l-^) (25) 
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The energy density p > follows a constraint G > —2FGf- The EoS of this YM field is 

where we have defined that 7 = e/G. When the condition of 7 = can be got at some state 
with i?^ 7^ and G{F) 7^ 0, the state of u; = — 1 is naturally realized. This condition can 
be easily satisfied. In the discussion as below we only consider these kind of YM fields. For 
example, in the model with the lagrangian ([3]), 7 = is got at the state E"^ = k,"^. Near this 
state, 7 > leads to u < —1, and 7 < leads to u > —1. But if the YM field has a trivial 
lagrangian with G = constant, which follows that 7 = 1, and u = 1/3. This is exactly the 
EoS of the relativistic matter, and it can not generate the state of 00 < 0. 

To study the evolution of EoS, we also consider the YM equation, which can be got by 
variation of S with respect to A^, 

d^{a^e F''^"') + r'^A^ia'^e F'^"") = 0, (27) 

from the definition of it is found that these equations become a simple relation: 

dr{a\E) = 0, (28) 

where E is defined by E"^ = Y^^^^Ef. If e = 0, this equation is an identity, and from ( l26|) . we 
know uj = —1, which can't be differentiated from cosmological constant. When e 7^ 0, this 
equation can be integrated to give 

a^eE = constant. (29) 

For we want to study whether or not the EoS of this YM field can cross u = —1, here we 
assume its initial state is ~ — 1. In this condition, from the expression of p and p, it 
follows that e ~ 0, and G{F) nearly keep constant, which is for the Universe is nearly 
de Sitter expansion and p ~ —G{F)E'^/2 is nearly a constant in this Universe. So the YM 
equation suggests that 

e oc a-^ (30) 



From the EoS of ( l26l) . one knows that 

a; + loca"l (31) 
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This is the EoS evolution equation of the general YM field dark energy models. It is exactly 
same with special case of Eq. fl2T|) . So it also keeps the characters of the special case with the 
lagrangian u will run to the critical point u = —1 with the expansion of the Universe. 
But it can not cross this critical point. These is the general characters of these kind of YM 
field dark energy models. For showing this more clear, we discuss two example models. 
First we consider the YM field with the running coupling constant 

G'(F) =5(F"-F,"), (32) 

where B and Fc are quantity with positive value, and n is a positive number. The constraint 
of p > follows that 

F > ^^=. (33) 
The dielectric constant can be easily get 

e = G + FG^ = 5(n + l)F"-5F,'^, (34) 

and 

7 = (n + l) + - V- (35) 



It is obvious that when F = F^/ \/n + 1, 7 = is satisfied, and which leads to 00 = —1. 
Near this critical state, E ~ F^/ (n + 1). So the YM equation of (l29l) becomes 



An 2„/ 2 
n + 1 \ n + 1 

which follows that 7 oc a^^. From the expression of uj in Eq. fl26l) . one can easily get 

u + l- -— oc (37) 

3 

This is exact same with the evolution behavior shown in formula flM]) . 

Another example, we consider the YM field with the coupling constant of 

G(F) = l-exp(l-|^), (38) 

where the constant quantity F^ 7^ 0. When F ^ F^, this lagrangian becomes the trivial case 
with G{F) = 1, but when F is near Fc, the nonlinear effect is obvious. Then 

e=l + (|^-l)exp(l-|^). 

-'^ c c 

so the critical state of F = 0.433Fc leads to 7 = and uj = —1. By the similar discussion, 
from the YM field (1291) . one can also get 7 oc near this critical state, which generates 
a; + 1 oc a~^. 
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III. TWO YM FIELDS MODEL 



In the former section, we have discussed that the dark energy models with single YM field 
can't form a state of uj crossing —1, not matter what kind of lagrangian or initial condition. 
But we should notice another character: the YM field has the EoS of oc — 1 + a~^, when 
its initial value is near the critical state of = — 1. So if the YM field has an initial state 
of > —1, it will keep this state with the evolution of the Universe, which is like the 
quintessence models. But if its initial state is < —1, it will also keep it, which is like the 
phantom models. This makes that we can build a model with two different free YM fields, 
one having an initial state of c<j > — 1 and the other being uj < —1. In this kind of models. 



111 ], where 



the behavior of uj crossing —1 is easily got. This idea is like the quintom models 
the authors built the model with a quintessence field and a phantom field. 

In the below discussion, we will build a toy example of this kind of model. Assume the 
dark energy is made of two YM fields with the effective lagrangian as Eq.Q 



L, = ^F,ln|^|, (2 = 1,2) (39) 
2 enf 



where Fi = E^{i = 1,2), and ki ^ K2- Their dielectric constants are 



From the YM field kinetic equations, we also can get the relations: 

ah^E, = Ci, (41) 

where Ci{i = 1, 2) are the integral constant, which are determined by the initial state of the 
YM fields. If the YM field is a phantom like field with Ui < —1, then < and Ci < 0. 
At the same time, a quintessence like YM field follows that Cj > 0. Here we choose the YM 
field of Li as the phantom like field with Ci < 0, and L2 as the quintessence like field with 
C2 > 0. The energy density and pressure are 

P^ = ^{e^ + b), p, = S(|_6), (42) 



SO the total EoS is: 



UJ = 



P1+P2 _ EKf-l) + El{f-l) 
p, + P2 El{(3r + l) + El{(32 + iy 



(43) 
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where we have also defined that (3i = ei/h. Using the relation of (3i and Ei, we can simplify 
the equation of state as 

3eft(/3i + l)a + e/52(/?2 + l)' ^ ' 

where a = ^1/^1- We need this dark energy has the initial state of > —1, which requires 
that the field of p2 is dominant at the initial time. This is easily obtained as long as at this 
time Ef{(3i + 1) < -E'|(/?2 + 1) is satisfied. The finial state, we need to get u < —1, which 
means that pi is dominant, and the behavior of crossing —1 realized in the intermediate 
time. But how to get this? From the before discussion, we know in the Universe with only 
one kind of YM field (i=l or 2), the YM equation follows that oc a'^. And it will go 
to the critical state ej = with the expansion of the Universe. At this state, Ef = k^, 
and Pi = hEf/2 = hn'j/2 keeps constant. So in this two YM fields model, if we choose the 
condition n\ > n\ {a > 1), this may follow that the finial energy density pi > p2, and pi 
is dominant. For this intent, we build this model with the condition as below: choosing 
a = 1.5, which can ensure the finial state, the first kind of YM field {i = 1) is the dominant 
matter. At the present time, corresponding to the scale factor ao = 1, we choose that 
Pi = —0.4 < and /?2 = 0.2 > 0, which keeps that the first field always having a state 
of cui < —1 (like the phantom) and the second field with > —I (like the quintessence). 
These choice of Pi leads to the present EoS 

4 e^i/3ia + 6^2/5 

^^=-1 + ^ a r = -1.10<-l, 45) 

which is like the phantom field. Since pi increases, and p2 decreases with the expansion of 
the Universe, there must exist of a time, before which p2 is dominant, and this leads to the 
total EoS > — 1 at that time. 

Combining the Eqs. fl41l) and (jBD, we can solve the evolution of EoS uo with the scale time 
in numerical calculation, where the relation of Ci and C2 is easily got 

For each kind of YM field, its EoS is 

"-ItI ^''^ 

The condition of > {Pi < 0) will generate uji > —1 {uji < —1). The evolution of them 
are shown in Fig. [3]. This is exact result we expect. Pi runs to the critical point of Pi = 

12 



with the expansion of the Universe, which makes Ui runs to a;^ = — 1, no matter what kind 
of initial values. This is same with the single YM field model. With Pi 0, the strength 
of the field Ef will also go to its critical point Ef — k}. This can be shown in Fig. [4]. For 
we have chosen the condition oi a = k\/ k\ > must lead to El > E^ at some time, the 
first kind of YM field becomes dominant, and the total EoS a; < — 1 is realized. This can 
be shown in Fig.[l]. In Fig. [2], we also plot the CMB power spectrum in the Universe with 
this kind of YM field dark energy, and find that it is difficult to be distinguished from the 
ACDM Universe, which is for the effect of the dark energy on the CMB power spectrum is 
an integral effect from CMB decoupling time to now. But the evolution detail of u is not 
obvious. This is the disadvantage of this way to detect dark energy. 

Now, let's conclude this dark energy model, which is made of two YM fields. One has the 
EoS of a;i < —1 and the other is a;2 > — 1. At the initial time, we choose their condition to 
make pi < P2, and second kind of YM field is dominant, which makes the total EoS u > —1 
at this time. This is like the quintessence model. For the ui < —1 keeps for all time, from 
the Friedmann equations, one knows that its energy density will enhance with the expansion 
of the Universe. And at last it will run to its critical point pi = hK\/2. And the same time, 
p2 will decrease to its critical state p2 = hK\/2. For we have chosen kI/k^ > which must 
make pi — p2 and some time, and after this, pi is dominant, the total EoS lu < —1. So the 
equation of state cu crossing —1 is realized. It is simply found that this kind of crossing must 
be from a;>— lto<— 1, which is exactly same with the observations. But the contrary 
condition, from a;<— lto>— 1 can't be realized in this kind of models. 

IV. CONCLUSION AND DISCUSSION 

In summary, in this letter we have studied the possibility of u crossing —1 in the YM 
field dark energy models, and found that the single YM field models can not realize, no 
matter what kind of their effective lagrangian, although this kind of models can naturally 
give a state of a; > —1 or a; < —1, which depends on their initial state. Near the critical 
state of a; — —1, the evolution of their EoS with the expansion of the Universe is same, 
0; + 1 oc a^^, which means that the Universe will be a nearly dc Sitter expansion. This is the 
most attractive character of this kind of models, and this makes it very like the cosmological 
constant. So the Big Rip is naturally avoided in this model. But this evolution behavior 
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also shows that the single field models can not realize uj crossing —1. This is same with the 
single scalar field models. 

But in these models, > — 1 and < — 1 all can be easily got. The former behavior is 
like a quintessence field, and the later is like a phantom field. So one can build a model with 
two YM fields, and one field with lo < —\ and the other with uj > —\. This idea is very 
like the quintom models. Then we give an example model and find that in this model, the 
property of crossing the cosmological constant boundary can be naturally realized, and we 
also found that this crossing must be from u;>— lto<— 1, which is exact the observation 
result. In this model, the state will also go to the critical state of u; = — 1 with the expansion 
of the Universe, as the single YM field models. This is the main character of the YM field 
dark energy models, which makes the Big Rip is avoided. The present models we discuss 
in this paper are in the almost standard framework of physics, e.g. in general relativity in 
4-dimension. There does not exist phantom or higher derivative term in the model, which 
will lead to theory problems in field theory. Instead, the YM field as ([3]), is introduced, 
which includes the gauge invariance, the Lorentz invariance, the correct trace anormaly, and 
the asymptotic freedom. These are the advantages of this kind of dark energy models. But 
these models also exist some disadvantages: first, what is the origin of the YM field? and 
why its renormalization scale is so low as the present density of the dark energy? In 
the two YM fields model, we must choose a > 1 to realize the lo crossing —1, which is a 
mild fine-tuning problem. All these make this kind of models being unnatural. These are 
the universal problems which exist in most dark energy models. If considering the possible 
interaction between the YM field and other matter, especially the dark matter, which may 



have some new character 
energy models 



27|. This topic had been deeply discussed in the scalar field dark 



but had not been considered in this paper. 
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FIG. 1: The evolution of the state equation w with the scale factor a. The model 1 denotes 
the single YM field model with the present EoS uiq = —1.2, model 2 denotes the model with 
(jOq = —0.8, and model 3 denotes the two YM fields model, where we have set kI = 1.5«;|, and 
Pi = —0.4, P2 = 0.2 at present with oq = 1. 
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FIG. 2: The CMB anisotropy power spectrum. The model 1 denotes the single YM field model 
with the present EoS ujq = —1.2, model 2 denotes the model with ojq = —0.8, and model 3 denotes 
the two YM fields model, where we have set = 1.5^2, and /3i = —0.4, /?2 = 0.2 at present with 
ao = 1. In this figure, the dots denote observation results from the WMApj^] satellite. Here we 



have used the CMBFAST program 
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FIG. 3: The evolution of the /?i and P2 with the scale factor a, where we have set /tf = I.SkI, and 
Pi = —0.4, P2 = 0.2 at present time with oq = 1. 




FIG. 4: The evolution of the Ei and E2 with the scale factor a, where we have set = 1.5k|, and 
Pi = —0.4, P2 = 0.2 at present with oq = 1. 
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